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Abstract 

Coherent structures are solutions to reaction-diffusion systems tliat are time-periodic in an appro- 
priate moving frame and spatially asymptotic at a; = ±oo to spatially periodic travelling waves. This 
paper is concerned with sources which are coherent structures for which the group velocities in the 
far field point away from the core. Sources actively select wave numbers and therefore often organize 
the overall dynamics in a spatially extended system. Determining their nonlinear stability properties 
is challenging as localized perturbations may lead to a non-localized response even on the linear level 
due to the outward transport. Using a modified Burgers equation as a model problem that captures 
some of the essential features of coherent structures, we show how this phenomenon can be analysed 
and nonlinear stability be established in this simpler context. 



1 Introduction 

In this paper, we analyse the long-time dynamics of solutions to the Burgers- type equation 

(j)t + ctanh(^Y) = ^xx + (pl, c>0 (1.1) 

with small localized initial data, where x G M, t > 0, and (j){x, t) is a scalar function. The key feature 
of this equation as opposed to the usual Burgers equation is that the characteristic speeds are c > at 
spatial infinity and — c < at spatial minus infinity: hence, transport is always directed away from the 
shock interface at x = and not towards x = as would be the case for the Lax shocks of the standard 
Burgers equation. 



We are interested in (1.1) due to its close connection with the dynamics of coherent structures that arise 



in reaction-diffusion systems 

ut = Du^^ + f{u), xEM, ueR"". (1.2) 
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Figme 1: Panel (i) shows the graph of a source u*{x,t) as a function of x for fixed time t: the group velocities 
of the asymptotic wave trains point away from the core of the coherent structure. Panel (ii) illustrates the 
Floquet spectrum of a spectrally stable source: the two eigenvalues at the origin correspond to translations 
in space and time, which, in contrast to the essential spectrum, cannot he moved by exponential weights. 
Panel (iii) shows the behaviour of small phase (p or wave number (px perturbations of a wave train: to leading 
order, they are transported with speed given by the group velocity Cg without changing their shape JT^. 



A coherent structure (or defect) is a solution u* {x, t) of ( 1.2 ) that is time-periodic in an appropriate moving 
frame y = x — c*t and spatially asymptotic to wave-train solutions, which are spatially periodic travelling 
waves of (1.2). Such structures have been observed in many experiments and in various reaction-diffusion 
models, and we refer to [6j for references and to Figure [T] for an illustration of typical defect profiles. For 
the sake of simplicity, we shall assume from now on that the speed c* of the defect we are interested in 
vanishes, so that the coherent structure is time-periodic. Coherent structures can be classified into several 
distinct types [21 |5l |6] that have different stability and multiplicity properties. This classification involves 
the group velocities of the asymptotic wave trains, and we therefore briefly review their definition and 
features. Wave trains of (1.2) are solutions of the form it (x, t) = u„t{kx — Lot; k) , where the profile ttwt(y; ^) 
is 27r-periodic in the y- variable. Thus, k and uj represent the spatial wave number and the temporal 
frequency, respectively, of the wave train. Wave trains typically exist as one-parameter families, where the 
frequency uj = u}^i{k) is a function, the so-called nonlinear dispersion relation, of the wave number k, which 
varies in an open interval. The group velocity Cg of the wave train with wave number k is defined as 

The group velocity is important as it is the speed with which small localized perturbations of a wave train 
propagate as functions of time t, and we refer to Figure [T|iii) for an illustration and to [Ij for a rigorous 
justification of this statement. The classification of coherent structures mentioned above is based on the 
group velocities Cg of the asymptotic wave trains at x = ibc«. We are interested in sources for which 
Cg < < c+ as illustrated in Figure [T|i) so that perturbations are transported away from the defect core 
towards infinity. Sources are important as they actively select wave numbers in oscillatory media; examples 
of sources are the Nozaki-Bekki holes of the complex Ginzburg-Landau equation. 

From now on, we focus on a given source and discuss its stability properties with respect to the reaction- 
diffusion system (1.2). Spectral stability of a source can be investigated through the Floquet spectrum of 
the period map of the linearization of (1.2) about the time-periodic source. Spectral stability of sources 
was investigated in [6] , and we now summarize their findings. The Floquet spectrum of a spectrally stable 
source will look as indicated in Figure [T|ii) . A source u*{x,t) has two eigenvalues at the origin with 
eigenfunctions M*.(x,t) and Uf{x,t); the associated adjoint eigenfunctions are necessarily exponentially 
localized, so that the source has a well defined spatial position and temporal phase. There will also be two 
curves of essential spectrum that touch the origin and correspond to phase and wave number modulations 
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Figure 2: The left panel contains a sketch of the space-time diagram of a perturbed source. The defect core 
will adjust in response to an imposed perturbation, and the emitted wave trains, whose maxima are indicated 
by the lines that emerge from the defect core, will therefore exhibit phase fronts that travel with the group 
velocities of the asymptotic wave trains away from the core towards ±00. The right panel illustrates the 



profile of the anticipated phase function 4){x, t) defined in (1.3) 



of the two asymptotic wave trains. It turns out that the two eigenvalues at the origin cannot be removed 
by posing the hnearized problem in exponentially weighted function spaces; the essential spectrum, on the 
other hand, can be moved to the left by allowing functions to grow exponentially at infinity. 

The nonlinear stability of spectrally stable sources has not yet been established, and we now outline why 
this is a challenging problem. From a purely technical viewpoint, an obvious difficulty is related to the fact 
that there is no spectral gap between the essential spectrum and the imaginary axis. As discussed above, 
such a gap can be created by posing the linear problem on function spaces that contain exponentially 
growing functions, but the nonlinear terms will then not even be continuous. To see that these are not just 
technical obstacles, it is illuminating to discuss the anticipated dynamics near a source from an intuitive 
perspective. If a source is subjected to a localized perturbation, then one anticipated effect is that the 
defect core adjusts its position and its temporal phase in response. From its new position, the defect will 
continue to emit wave trains with the same selected wave number but there will now be a phase difference 
between the asymptotic wave trains at infinity and those newly emitted near the core. In other words, 
we expect to see two phase fronts that travel in opposite directions away from the core as illustrated in 
Figure [2| The resulting phase dynamics can be captured by writing the perturbed solution u{x,t) as 

u{x, t) = u*[x + (j){x, t), t) + w{x, t), (1-3) 

where we expect that the perturbation w{x,t) of the defect profile decays in time, while the phase (j){x,t) 
resembles an expanding plateau as indicated in Figure [2] whose height depends on the initial perturbation 
through the spatio-temporal displacement of the defect core. 

The preceding discussion indicates that localized perturbations of a defect will not stay localized but result 
instead in phase fronts that propagate towards infinity. As a first step towards a general nonlinear stability 
result for sources in reaction-diffusion systems, we focus here on the Burgers-type equation 



-|- c tanh 



+ 



/)2 







(1.4) 



for the phase function 4'{x,t) introduced in (1.3). The rationale for choosing this model problem is as 



follows. First, as established formally in [3j and proved rigorously in [T], the integrated viscous Burgers 
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equation captures the phase dynamics of wave trains over long time intervals. Since the Burgers equation 
does not admit sources, we add the artificial advection term on the left-hand side which creates the 
characteristic speeds itc at x = ±00 that account for the outgoing group velocities of the asymptotic 
wave trains of the source. While the inhomogeneous advection term can be thought of as fixing the 
position X = of the core, the equation still has a family of constant solutions, which correspond to 
different temporal phases of the underlying hypothetical sources. We therefore feel that gaining a detailed 



understanding of the long-time dynamics of (1.4) for small localized initial data will shed significant light 



on the expected dynamics of sources and on the techniques needed to analyse their stability. We emphasize 
that the dynamics of wave trains of reaction-diffusion systems under non-localized phase perturbations was 
investigated only recently in j7]; the methods used there rely on renormalization-group techniques which 
seem difficult to generalize to sources. On the other hand, the analysis presented here is currently far more 
limited in terms of the equations it applies to. 



The linearization of (1.4) about = is given by 



ctanh 



(1.5) 



The spectrum of the operator in (1.5) is as shown in Figure [Tj^ii) except that there is only one embedded 
eigenvalue at the origin that cannot be moved by exponential weights: the eigenfunction of the eigenvalue 
at A = is given by the constant function |, and the associated adjoint eigenfunction is given by 



tp{y) = sech^ ^ 



Equation (1.5) admit^the Green's function 

1 (x-y + ct)^ 



Q{x,y,t) 



+- 



errfn 



1 



+ 



{x-y-ct) 
At 



1 



y — X + ct 



it 



where the error function is given by 



errfn(z) 



1 



errfn 



ds. 



1 + e-'^y 

X — ct 



(1.6) 



it 



V'(y), 



The first two terms in the Green's function are Gaussians that move with speed itc away from the core 
and decay at the rate l/\/t, while the term comprised of the difference of the two error functions produces 
a plateau of constant height | that spreads outward as indicated in Figure [3[ 

For sufficiently localized initial data (j)o{x), the solution (p{x,t) to the linear equation (1.5) is therefore 
given by 

4>{x,t) = / g{x,y,t)(j)o{y)dy, 



which converges to a constant: 



(/)(x,t) 



V'(y)0o(y)dy as t 



00 



^Note that ( 1.5 1 is the formal adjoint of the hnearization of the standard viscous Burgers equation ut — it. 
(c/2)[l - tanh(c2;/2)] with x = 



2uuz about 

ct, whose Green's function can be found via the hnearized Cole- 



the Lax shock u{x) 

Hopf transformation by setting w{x, t) = cosh(cx/2) u{y, t) Ay. The Green's function of ( 1.5 1 can then be constructed by 
reversing the roles of x and y in the Green's function for the Lax shock linearization. 
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Figure 3: Shown are the graphs of the function errfn((— z + ct) /-s/it) — errfn((— 2; — ct) / -\/4t) for smaller and 
larger values of t, which resemble plateaus of height approximately equal to one that spread outwards with 
speed ±c. 



for each fixed x. In particular, one cannot expect the solution associated with localized initial data to 



remain localized for all time. In fact, the same is true for the nonlinear equation (1.4): the Cole-Hopf 
transformation 



{x,t) _ 



{x, t) = log 



l + (t){x,t) 



relates solutions (j){x,t) of (1.4) and solutions 4>{x,t) of the linearization (1.5), and we conclude that 



solutions (j){x,t) of (1.4) are given by 



1 + / Gix,y,t)(j)o{y)dy 



(x, t) = log 



For t — )• 00, these solutions converge again pointwise in x to the constant 



log 



The Cole-Hopf transformation does not, however, extend to more general equations, and we therefore 
pursue here a different approach that, we hope, will also be useful when investigating the nonlinear stability 
of sources in general reaction-diffusion systems. 



In order to analyze the dynamics of (1.4) in a way that may also be applicable in the case of a general 



reaction-diffusion equation, we need to find an appropriate ansatz for small-amplitude solutions. To do so, 
we define 13{x, t) by 

B{x,t) := g(x,0,t + 1), 



where Q{x,y,t) is the Green's function defined in (1.6), and note that 

qy*{x,t,p):=log{l+pB{x,t)) (1.7) 
is then a solutiorj^of (1.4) for each fixed p G M. We now investigate the long-time dynamics of solutions 



to (1.4) with initial data (l){x,0) = (t)o{x) using the ansatz 

^{x,t)=(P*{x,t,p{t))+v{x,t), (1.8) 
where p{t) is a real-valued function, and v{x, t) is a remainder term. At time t = 0, we normalize the 



decomposition in (1.8) by choosing p{0) = po such that 



ip{x)[(po{x) - (j)*{x,0,po)] dx = 0. 



In fact, (j)*(x,t,p) does not need to satisfy (1.4 1 exactly: our analysis goes through provided B{x,t) can be chosen such 



that (ji*{x, t,p) satisfies ( 1.4| approximately with an error of order of (1 + ^ times a Gaussian. 



5 



We will prove in ^ that a unique pq with this property exists for each sufficiently small localized initial 
condition (pQ. Note that the leading behavior (j)*{x, t,p{t)) corresponds to a plateau with height log(l+p(t)) 
of length 2ct that spreads outward with speed itc. The main result of this paper is as follows. 

Theorem 1. For each 7 G (0, there exist constants eo,rio,Co, Mq > such that the following is true. 
If 4>Q G C"^ satisfies 

e:=||e"'/^°(/>o||ci <eo (1.9) 



then the solution 4>{x, t) of (I.4) with (f){-, 0) = (f>o exists globally in time and can be written in the form 

(f){x, t) = (f>*{x, t,p{t)) + v{x, t) 



for appropriate functions p{t) and v{x,t) with (p* as in (1.1). Furthermore, there is a poo G M with 
\Poo\ < Co such that 

W) - Pool < eCoe-^°\ 

andv{x^t) satisfies 

\v{x,t)\ < j^-^ (^e"Wi) +e"s?o(^+i)j , \v,{x,t)\ < j—^^{e~^'^+e-^Wi+^ 



for all f > 0. In particular, \\v{-,t)\\Lr — as t — )• 00 for each fixed r > 

Note that our result assumes that the initial condition is strongly localized in space. We believe that this 
assumption can be relaxed significantly; for the purposes of this paper, however, phase fronts are created 
even by highly localized initial data, and the key difficulties are therefore present already in the more 
specialized situation of Theorem [TJ 

The exponential convergence of p{t) reflects the intuition that a source has a well-de&ied position due to 
the exponential localization of the adjoint eigenfunction ip. The asymptotics of the perturbation v is given 
by moving Gaussians that decay only like (1 + 1)^'^ for each fixed 7 G (0, ^), rather than with (1 + 1)~ 2 as 



expected from the dynamics of the viscous Burgers equation. This weaker result is due to the form (1.8) 
of our ansatz, which effectively creates a nonlinear term that is proportional to g{x)uux for some function 
g{x). While this term resembles the term 2uux in Burgers equation, it does not respect the conservation- 
law structure. Thus, we only obtain decay at the above rate. Although this may not be optimal, it allows 
us to avoid terms that grow logarithmically in the nonlinear iteration in ^ We do not know if it is possible 
to improve this rate to 7 = ^ by adjusting our ansatz appropriately. 

The remainder of this paper, ^ is devoted to the proof of Theorem [T] 



2 Proof of the main theorem 

To prove Theorem [l| we set up integral equations for {p, v) and solve them using a nonlinear iteration 
scheme. Throughout the proof, we denote by C possibly different positive constants that depend only on 
the underlying equation but not on the initial data or on space or time. 
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2.1 Derivation of an integral formulation 



Substituting the ansatz 

t) = log (1 + p{t)B{x, t)) + v{x, t), B{x, t) := g{x, 0,t + l) 
into the Burgers-type equation 



(pt + c tanh 

we find that {p, v) needs to satisfy the equation 



cx 



(CX \ p 
— j Vx + v'^- ^ g(x,0,t + 1) +Af{x,t,p,p,Vx), 



(2.1) 



where the nonhnear function M is given by 



.r, , . ^ 2pvA{x,t) . f B{x,t) 
M{x,t,p,p,Vr,) := , , +p 



B{x,t) 



l+pB{x,t) '"'yi + lp l+pB{x,t))' 
The idea of the proof is to use an appropriate integral representation for v that will allow us to set up a 



nonlinear iteration argument to show that solutions (p, f) of (2.1) exist and that they satisfy the desired 
space anc 



decay estimates in space and time. Recall from (1.6) the expression 

.2 I 



1 (x-y + cty 

e « 



+ - 



/Ivrt 
c 



^errfn ^ 



1 + e'^y 
y — x + ct 



1 _ (x-y — ct) 

H ^=e 4i 



1 



4t 



Ant 1 + e-^y 

y — X — ct 



(2.2) 



errfn 



4t 



for the Green's function of the linear problem (1.5) and note that it satisfies the identity 



/ 



g{x, y, t - s)g{y, 0, s + 1) dy = g{x, 0, t + 1). 
Using this identity and the variation-of-constants formula, we can rewrite (2.1 ) in integral form and obtairj^ 



v{x, t) 



cp{0) 



^g(:r,0,t + l)(log(l + ^)-log^l+ ^ 

t 



)) 



(2.3) 



+ / 0{x,y,t)vo{y)dy + / / g{x,y,t - s) [vy + ■,P,P,Vy)] {y, s) dy ds. 



The expression (2.2 ) of the Green's function g{x, y, t) shows that the terms involving the error functions do 



not provide any temporal decay. We need to treat these terms separately and therefore write the Green's 
function g{x,y,t) as 

0{x,y,t) = £{x,y,t) + g{x,y,t), 



where 



Six, y, t) = e{x, t)'ip{y), e{x, t) := 



^errfn ^ 



X + ct 



errfn 



X — ct 



^We shall often use the notation [/ + g]{x,t) to denote f(x,t) + g{x,t). 
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and 



Gix,y,t) := g{x,y,t) - £{x,y,t) 



(2.4) 



1 (x^y + ct)^ 1 1 jx-y-ct)'^ 
e 4t z — ^ — H ; P 4t 



/47rt 



1 + e^?' Vivrt 



1 + e-'^y 



+ ( errfn f ^ ^ | — errfn 



4t yy 4 ^2' 



/ ^ (-X + Ct\ f-X-ct\\C 

errm — — — errm ' ' ' " — ' 



4t 



At 



^sech { — ). 



A calculation shows that 
' y — X ± ct' 



errfn 



4t 



errfn 



/— X ± ct 



V V4t 



I sech2(y) < Ct"^/^ f e"'" 4!"' + e « |e 



■c|y|/4 



and we conclude that 



\G{x, y,t)\< Ct-'/^ e--^^ + e 



(x — y-ct) 
At 



(2.5) 



To analyse (2.3), we consider the initial condition (j)o{x) = 0*(x,O,p(O)) + vq{x) and show that, if cfiQ is 



sufficiently small in L°°, then p(0) = po can be chosen such that 



i'{y)vo{y) dy = 



or, equivalently. 



i^{y)[My)-<P*{y,o,po)]dy = o. 

To prove this claim, we observe that 

(/<*(2/,0,po) = log(l +Po^(2/,0)) =poa(y,0, 1) + 0(p2). 



(2.6) 



(2.7) 



We note that the term 0(pq) in (2.7) is bounded uniformly in y G M, and substitution into (2.6) therefore 
gives the equation 

V'(y)'Ao(y)dy = po / i^{y)G{y,o,i)dy + 0{pl), 



which can be solved uniquely for pQ = p(0) near zero for each (pQ E L°° for which ||(/'o||l°° is small enough. 
In particular, there is a constant C > such that the resulting initial value p{0) satisfies 

|p(0)| < C||0o||l- <Ce. 

We will construct p{t) such that 

Pit) = (1 + ^) ^V'(y) [vl+M{;;P,P,Vy)] {y,t)dy (2.8) 

or, equivalently, that 

log (1 + ^) = log (1 + ^) + 1 1^ I^Hy) H +AA(-, ■,p,p,vy)] {y,s)dyds (2.9) 
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for all t>0. Substituting (2.9) into (2.3), we obtain the equation 
v{x,t) = 



-g{x, 0,t + l) (log (l + ^) - log f 1 + ^) ) + / g{x, y, t)vo{y) dy (2.10) 



4 7 



4 J 



+ / / G{x,y,t- s)[vl+M{-,-,P,P,Vy)]{y,s)dyds 

Jo Jm. 

+ ie{x,t-s)-e{x,t + l))^p{y)[vy+^^i■,■,P,P,Vy)]{y,s)dyd^ 

Jo Jm 



for v{x, t). 



To analyse solutions of (2.8) and (2.10), we introduce template functions that capture the anticipated 
spatio-temporal behaviour of solutions. For each fixed choice of 7 G (0, 5) and M > 0, we define 



i{x,t) 



{i + ty 



g M{t+1) _|_ g M(t+1) 



ix,t) 



(l + t)7+l/2 



g M(t+1) _|_ g M{t+1) 



(2.11) 



and let 



hiit) :-- 



sup 

,0<s<t 



h 02 



iy,s), 



h2{t) := sup |p(s)|e^'^/*^ h{t) := hi{t) + h2{t). 

0<s<t 



We will later choose M » 1. 

We remark that we do know existence and smoothness of {v,p) for short times: Indeed, we can solve the 



original PDE for (j)[x,t) for short times and can substitute the resulting expression into (2.8) upon using 



(1.8). The resulting integral equation has a solution pit) for small times and, using again (1.8), we find a 



smooth function v that then satisfies (2.10). Furthermore, using that satisfies (1.9) by assumption, we 



see that h{t) is well defined and continuous for < t <C 1. Finally, standard parabolic theory implies that 
h{t) retains these properties as long as h{t) stays bounded. The key issue is therefore to show that h{t) 
stays bounded for all times t > 0, and this is what the following proposition asserts. 



Proposition 2.1. For each 7 E (0, ^), there exist positive constants eo,Co,M such that 

hi{t) < Co(e + h2{t) + h{tf), h2{t) < Co(e + h{tf). 
for all t > and all initial data uq with e := \\e^^^^^ (f)Q\\(ji < eo- 



(2.12) 



Using this proposition, we can add the inequalities in (2.12) and eliminate /12 on the right-hand side to 
obtain 

h{t) < Co{Co + l){e + h{tf). 

Using this inequality and the continuity of h(t), we find that h{t) < 2Co(Co + l)e for all t > provided 
< e < eo is sufficiently small. Thus, Theorem [T] will be proved once we establish Proposition 2.1 The 
following sections will be devoted to proving this proposition. 



2.2 Estimates of the nonlinear term 



We begin by deriving estimates of the nonlinear term 



. . 2pv^Br,{x,t) ./B{x,t) B{x,t) 

Al{x,t,p,p,v^) = — — — — — +p ' 



l+pB{x,t) "^Vl + fP l+pB{x,t) 
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that appears in (2.8) and (2.10). Recalling that 

e(x, t) = T I errfn ( , 

4 V V \/4t 

and using the definition of the error function, we see that 



X + ct\ f X — ct 

errfn 

At 



eix,t) ( - - e{x,t) 



< C e 8t + e 8t 



Since B{x, t) = e{x, t + 1) + G{x, 0,t + 1), we conclude from (2.5) that 



and thus 



In addition, we have 



B{x,t)(--B{x,t) 



B{x,t) B{x,t) 



< C e + e 



l+pB{x,t) 1 + Ip 



< C\p\ e + e 8(t+i) 



B^{x,t) 



1 +pB{x,t) 

Combining these estimates, we therefore obtain 



< C(l + ty^/"^ ( e + e 



\Af{x,t,p,p,v^)\ < C ({l + ty^/^\p\\v^\ + \pp\) ( e"^(W +e-^(wT 



uniformly in a; E M and t >0, provided p is sufficiently small. 



(2.13) 



2.3 Estimates for h2(t) 



To establish the claimed estimate for h2{t), recall from (2.8) that 

\pit)\<Cii + \pm [ i,{y) [vl+M{-r,V,V.Vy)]{y,t) 

JR 

We shall show that there is a constant Ci = Ci{M) such that 

m\<Cie-'^"'/^{e + h\t)), 



dy. 



(2.14) 



(2.15) 



which then establishes the estimate for /i2(t) stated in Proposition 2.1 To show (2.15), we note that the 
definitions of p{t) and /i2(i) imply that 



m\ < 



\pm+ ['\p{s)\ds<\pm+ [\-'^''/^h2{t)ds<C,{e + h2{t)) (2.16) 
Jo Jo 



and therefore 



\p{t)pit)\ < C7ie-^'*/*^/i2(t)(e + /i2(t)) < Cie-'''/^'ie + hitf). 



(2.17) 



Next, we use the estimate |'0(y)| ^ 2e '^'^l, the bound (2.13) on J\f, and the inequality 



c\y\ («±ct) 



e 2 e M(i+t) < (7^g 4 Q M 
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which holds for each M > 8, to conclude that 



(2.18) 



and therefore 



\^{y)Af{y,t,p,p,Vy)\ < Ci(l + t)-i/2 's^+i) + e ) V'(y)l%(y, OIIpWI 

Q-^(J+TT _^Q-^(t+TT \ 'il){y)\p{t)p{t)\ 

< Ci{l + t)~^/^e-l\y\-^'hi{t){e + h2{t)) + Cie-f ls'l-Tr*(e + h{tf) 

< Cie-^2\y\~^\e + h{tf). (2.19) 



Using these estimates in (2.14), we arrive readily at (2.15), thus proving the estimate for /i2(t) stated in 



Proposition 2.1 



2.4 Estimates for v{x,t) and Vx{x,t) 
In this section, we will estabhsh the pointwise bounds 

\vix,t)\ < Ci{e + h2{t) + h{tf)0i{x,t) 
\vx{x,t)\ < Ci{e + h2{t) + h{tf)e2{x,t) 

for v{x,t) and Vx{x,t), respectively, which taken together prove the inequality for hi{t) stated in Proposl 



(2.20) 
(2.21) 



tion|2.1[ In particular, the proof of Proposition 2.1 is complete once the two estimates above are established. 



We denote by Ci possibly different constants that depend only on the choice of 7, M so that Ci = Ci(7, M). 



We focus first on v{x,t). From the integral formulation (2.10) for v{x,t), we find that 
\vix,t)\ < Cii\po\ + \p{t)\)gix,0,t + l)+ [ g{x,y,t)\voiy)\dy 

dy ds 



(2.22) 



f I G{x,y,t- s)[vl+J\f{-,-,p,p,Vy)]{y,s) 

Jo JR 



+ 



z{x,t-s + l) - e(x,t + 1)]V'(2/) [vy +Af{-, Vy)] {y,s) dyds 



In the remainder of this section, we will estimate the right-hand side of (2.22) term by term. 



First, we note that (2.5) implies that there is a constant Ci = Ci{M) so that g{x,0,t + 1) < Ci9i{x,t). 



Using (2.16) and the fact that \p{0)\ < Cie, we therefore obtain 

(|p(0)| + \pmg{x, 0, t + 1) < C7i(e + h2{t))ei{x, t) 



which is the desired estimate for the first term on the right-hand side of (2.22). 



Next, we consider the integral term in (2.22) that involves the initial data vq. Using (2.5) together with 



our assumption ( |1.9 ) on (/>o, and hence on vq, we see that 
\gix,y,t)voiy)\dy<Ce [ r^l^ 



J 



(x-y + ct) (x-y-ct) \ y~ 

e 4t + e 4t I e Af dy. 



(2.23) 
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which is clearly bounded by eCi9i{x, t) for t > 1 upon using 



y \ — ^ L 1^ y y 



For f < 1, we can use the estimates 



-1/2 



8f e < 2e »t Q M < C\e 2m 



(3:-i; + ct)^ 



_ (x — y — ct) 



St -)- e 8t 



dy < Ci 



to conclude that the integral in (2.23) is again bounded by eCi9i{x,t) 



We now consider the remaining two integrals in (2.22). Note that the definition of hi gives 



\vy{y,s)\ < e2{y,s)hi{s). 



Using this fact together with the estimates (2.16) and ( |2.17 ) and the bound (2.13), we obtain 



\Ar{y,s,p{s),p{s),Vy{y,s))\ < Ci(l + s)"i/2 e 8(s+i) + e \vy{y, s)\\p{s)\ 



{y+csy 



+Ci(e sf'^+i) +e 8(^+1) \\p[s)p{s)\ 
< Ci(e + h2{t))hi{t)e2{y, s)(l + s'y^'^ (e'^-^W^ + e"^^!^ 



+Ci(e + /i(t)^)e"Ar e 8{.+i) + e 



< Ci{e + h{tf) 



(1 + s)7-i/2^,(y, s)e2{y, s) + {l + s)e-ivr 



In ^2.5, we will prove the following result for our spatio-temporal template functions. 
Lemma 2.2. For each sufficiently large M , there is a constant Ci so that 




\g{x, y,t-s)\ [ej + (1 + sy-^/^eie2 + (i + sroie-'''/^' 
mx,y,t- s)\ k + (1 + sr-'/'ei02 + (1 + syeie-^'^/^' 



{y,s)dyds < Ci6i{x,t) 
{y,s)dyds < Ci02{x,t). 



Using this lemma and the above estimates for Vy + A^, we obtain the desired estimate 

g{x,y,t-s) [vl+M{; ■,P,P,Vy)] {y,s)dyds < Ci{e + h{t)'^)6i{x,t). 
Finally, we have the following lemma, whose proof is again given in §2.5[ which provides the desired estimate 




for the last integral in (2.22). 



Lemma 2.3. For each sufficiently large M , there is a constant Ci so that 




^{x,t-s + l)-e{x,t + l)\ij{y) [vl+M{-,-,P,P,Vy)] iy,s)dyds < Ci(e + /i(t)2)^i(x, 
3^{x,t - s + 1) - e^{x,t + l)\tP{y) [v^+Af{-,-,p,p,Vy)] {y,s)dyds < Ci{e + h{tf)92{x,t). 
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In summary, combining the estimates obtained above, we have estabhshed the claimed estimate (2.20), and 



it remains to derive the estimate (2.21) to complete the proof of Proposition 2.1 Taking the x-derivative 



of equation (2.10), we see that 



\vAx,t)\ < i\p{0)\ + \p{t)\)g^{x,0,t)+ g^{x,y,t)voiy)dy 



+ 



(2.24) 



Qx{x,y,t- s) [vl+N{-,-,P,P,Vy)\ {y,s)dyds 

+ / [ex{x,t-s + l)-ea:{x,t + l)]ip{y)[vy+M{-,-,P,P,Vy)]{y,s)dyds. 
Jo Jr 

Applying the second estimate in Lemmas 



2.2 



and 



2.3 



to (2.24) and using that g^{x,0,t) < C92{x,t), we 



immediately obtain (2.21). 



2.5 Proofs of Lemmas [272] and Q 



It remains to prove the lemmas that we used in the preceding section. 



Proof of Lemma 2.2. We need to show that for each large M there is a constant Ci so that 



g{x,y,t-s)\ [^i + (1 + 5)^-1/2^1^2 + (1 + s)^^ie-'='"/*^] {y,s)dyds < Ci0i(x,t) 

for all t > 0. First, we note that there are constants Ci, Ci > such that 

Cie-y"/'' < \9i{y,s)\ + \e2{y,s)\ < Cie-^''/^ 

for all < s < t < 1. Thus, for some constant Ci that may change from line to line, we have 

ft 



\g{x,y,t- s) 



+ a + s 



>7-l/2/ 



h + a + sVOie-^ 





< Ci [ /"(t-s)-i/V^^e~M dyds 
Jo Jr 
ft 



{y,s)dy ds 



< Ci 







{x — y) JO 

(t - s)-i/2e" W^e"8(^ dy + / {t - sY'^^'^fT^^^^ dy 

{\y\>'i\A} Am<2N} 



ds 



e S(t-a) -|- g M 



ds 



< Cie M 

< ^Oi{x,t) 

for all < t < 1. An analogous computation can be carried out for the x-derivative since (t — s)~^^'^ ds 
is bounded uniformly in < t < 1. 

Thus, it remains to estimate the expression 

ft 



;i(x, t)-' I / \g{x, y,t-s)\ ei + (1 + sr-'^'ei92 + (i + sye^e 

10 Jr 



{y,s)dyds 



13 



for t > 1. Combining only the exponentials in this expression, we obtain terms that can be bounded by 

{x + asct)'^ {x — y + aic{t — s))'^ {y + a2cs)'^ 



exp 



(2.25) 



M{l + t) 4(t-s) M(l + s) 

with Uj = ibc. To estimate this expression, we proceed as in ^ Proof of Lemma 7] and complete the 
square of the last two exponents in (2.25). Written in a slightly more general form, we obtain 

{x - y - ai{t - s))'^ {y - a2s)'^ {x - ai{t - s) - a2s)'^ 



Mi{t - s) 



+ 



M2(l + s) Mi{t - s) + M2{1 + s) 



^ Mi{t - s) + M2{1 + s) _ xM2{l + s) - (aiM2(l + s) + a2Mis){t - s) 



MiM2{l + s){t - s) 



Mi{t - s) + M2{1 + s) 



and conclude that the exponent in (2.25) is of the form 

(x + aat)^ {x — ai{t — s) — a2s)'^ 
M{l + t) A{t- s)+M{l + s) 



(2.26) 



4(t - s) + M{l + s) 
AM{l + s){t- s) 



y 



xM{l + s)- (aiM(l + s) + 4a2s)(t - s] 
4(t-s)+M(l + s) 



with aj = ibc. Using that the maximum of the quadratic polynomial ax"^ + /3x + 7 is (4a) + 7, it is 
easy to see that the sum of the first two terms in (2.26), which involve only x and not y, is less than or 
equal to zero. Omitting this term, we therefore obtain the estimate 

(xibct)^ (x - y5ic(t - s))^ {y-b2csf^ 



exp 



< 



M(l + t) 4(t-s) 
4(t-s) + Ms 



exp 



4M(l + s)(t-s) 



M(l + s) 

xM(\ + s) + c((5iM(l + s) + 452s)(t - s) 
4(t - s) + M(l + s) 



(2.27) 



for bj = ±1. Using this result, we can now estimate the integral (2.25) term by term using the key 
assumption that < 7 < ^. The term involving can be estimated as follows using (2.27): 



(x,t) 



-1 



\Q{x,y,t- s)\el{y,s)dyds 



Jo 



X / exp 



< Ci(l + t)" 



^/t^(l + s)l+27 

4{t- s)+M{l + s) 
AM{1 + s){t- s) 



[xM{l + s)± c(M(l + s) + 4s) (t - s) 
4(t-s)+M(l + s) 



dy ds 



AM{l + s){t-s) 



Vt^{l + s)i+27 W 4{t -s) + M(l + s) 



ds 



t/2 



^ds + Ciil + t)' 



(l + s)l/2+27 (l + t)l/: 
< Ci(l + t)^-i/2^Ci(l + t)-^ 

which is clearly bounded since 7 < ^. Similarly, we have 



t/2 (1 + ^)1+27 



ds 



i{x,t) 



-1 



\g{x, y, t - s) I (1 + sy-'/^9ie2 {y, s) dy ds 
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< Ci{i + ty 



1 / 4M{l + s){t-s) 

Vt^il + sj^ V 4(t -s) + M(l + s) 

t/2 1 



ds 



< C7i(l + t)^-i/2 + Ci, 



which is again bounded due to 7 < i. Finally, we estimate 



lix, ty 



< Ci(l + t)" 



\g{x,y,t- s)\{l + s^e M 6*1 (y, s) dy ds 



e M 



4M{l + s){t- s) 



y 4(t - s) + M(l + s) 

*/2 2. 



ds 



< Ci{l + tr I e-V \ ds + Ci{l + tre-^ / ds 

< Ci(l + t)^-i/2 + Ci(l + t) 



e 2M J 



which is bounded, again due to 7 < I. 



It remains to verify the second inequality in Lemma 2.2 which involves Gx- We shall check only the term 
involving (1 + s)"'~^^'^6i92 as the other cases are similar and, in fact, easier. We have shown above that 
the resulting integrals are bounded for < t < 1 and therefore focus on the case t > 1. Using that 



\Gx\ < Ct ^^'^\G\, which follows by inspection, and employing again (2.27), we obtain 
e2{x,t)-^ f [ \G.{x,y,t-s)\{l + sy-^/^Bi92{y,s)dyds 

Jo JR 



Jo 
ft 

< C7i(l + t)T+i/2 / 

Jo 

< Ci{i + tyt-^/^ + c\, 



4M(l + s)(i-s) 



(t - s)(l + s)^+i V 4(t - s) + M(l + s) 



ds 



tV2(l + s)7+l/2 (l+t)l/2 



ds + Ci 



t/2 (t-s)l/2(l + t)l/2 



ds 



which is bounded for t > 1. This completes the proof of Lemma 2.2 



□ 



Proof of Lemma 2.3 We need to show that 




\e{x,t- s + 1) - e{x,t + + AA(-, ■,p,p, Vy)]{y, s) dy ds < C7i(e + h{tf)ei{x,t). 



Intuitively, this integral should be small for the following reason. The difference e{x,t — s) — e{x,t + 1) 
converges to zero as long as s is not too large, say on the interval s G [0,t/2]. For s G [t/2,t], on the other 
hand, we will get exponential decay in s from the localization of i/jiu) in combination with the propagating 
Gaussians that appear in the nonlinearity and forcing terms. To make this precise, we write 



e{x, t — s) — e{x, t + 1) = e{x, t — s) — e{x, t — s + 1) + e{x, t — s + 1) — e{x, t + 1) . 



(2.28) 



term I 



term II 
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We focus first on the term I and consider the cases t > 1 and < t < 1 separately. First, let t > 1. For 
< s < i — 1, we have \e{x, t — s) — e{x, t — s + 1)\ < CQ{x, 0,t — s), and we can estimate the resulting 
integral above in the same way as in the proof of Lemma |2.2[ we omit the details. For t — l<s<t<l, 
on the other hand, the definition of e{x, t — s) yields 



\e{x, t — s) — e{x, t — s + 1)\ < C 



(l + t-s) 



Using (2.18)-(2.19), namely 



rnvWiy, s,p,p, vy)\ < Ce-l\y\-'-t^{e + h{tf) 



(2.29) 



(2.30) 



for all s > 0, we obtain 



[e{x, t- s) - e{x, t- s + l)]ip {y)\J\f{y, s,p,p, Vy)\{y, s) dy ds 

t-i 

2 o 2 2 „ 2, 

< Ci{e + h{tf) / e-^e-Tds < Ci{e + h{tf)e-^ e'^ , 



which is clearly bounded by Ci6i{x,t) since e < Ci(l + 1) ^ and 

(x + ct)'^ 2x^ 4c^ ch 
< — + ^ + 



M{l + t) - M M M 



for arbitrary M > 4. In summary, we have established the desired estimates for the term I in (2.28) for 



t > 1. For t < 1, the estimate (2.29) remains true since t — s is small, and proceeding as above yields 
ft 

[e{x,t- s) - e{x,t- s + l)]ij{y)\M{y,s,p,p,Vy)\dyds < C {e + h{tf)e~ e~ , 







which is again bounded by Ci0i(x,t). 

It remains to discuss the term II which involves the difference e{x^t — s + 1) — e{x,t + 1). We have 



|e(a;, t — s + 1) — e{x, t + 1 
rt-s+l 

er{x, r) dr 

't+i 
ft+i 



< 



t-s+l 

t+1 



V4 



vrr 



e 4t -)- e 4t 



+ 



{x — Ct) (x-ct)^ (^x + CT) (x+CTf 

e 4t — e 4t 



4r 



4t 



dr 



lt-s+1 \ Vt T 



{x — ct) 



+ e" 



(x + ct)^ 



dr, 



where we used in the last inequality that ze ^ is uniformly bounded in z. We now use the preceding 
expression to estimate ^j~^(x, t){e{x, t — s + 1) — e{x, t + 1)) and focus first on the single exponential term 



(x-ct)^ (x-ct)^ 
g M(l + t) g St 



Combining these exponentials and completing the square in x in the resulting exponent, the latter becomes 



[M{t-T + 1) + {M -8)t] 
8M(t + l)r 



x + 



c(8 - M)r(t + 1) 
M{t - r) + (M - 8)r 



1 2 



+ 



C^t-T + l) 



M{t-T + 1) + {M -8)t' 
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Using that t < t and picking M > 8, we can neglect the exponent resulting from the first expression that 
involves in x and conclude that 



eM(i+t)e 8^ < Cie" 



The remaining exponentials can be estimated similarly, and we obtain 



l^\x,t)\e{x,t-s + l)-e{x,t + l)\ < Ci{l + ty['^ ( ^ + -) 

Jt-s+i Wt tJ 



e M dr 



< Ci{l + ty{l + t-s)~^/^ei^. 



Using this inequality together with (2.30) finally gives 

0i{x,ty^ I I [e{x,t - s + 1) - e{x,t + l)]^l){y)\J\f{y,s,p,p,Vy)\{y,s)dyds 
Jo Jr 



2 rt 



< Ci{l + t)^{e + h{ty) / (l + t-s)-^/^e^e-^ds 

r 2 

< Ciil + ty{e + h{tf) (l + t)-i/2/ e-%/ds + e-fA7 / {1 + t - s)-^^"^ di 

Jo Jt/2 

< Ci{e + h{tf). 



for M sufficiently large, which proves the first estimate in Lemma 2.3 
It remains to prove the estimate 

- s + 1) - ea,{x,t + l)\'il){y) [vl + J\f{-,-,p,p,Vy)\ {y,s)dyds 
< Ci{e + h{tf )62{x,t) 

for the derivative in x. Since the derivative ofe(x,t — s + — e{x^t + 1) with respect to x generates an 
extra decay term (1 + t)"-'^/^, we have 

6'2(x,t)~i / / [ex{x,t- s + l) - e^{x,t + l)]tlj{y)\N{y,s,p,p,Vy)\{y,s)dyds 
Jo Jr 

< Ci(e2 + /i(t)2)(l + t)7+V2 /■ (i + i_s)-ie§^e-^ds 

Jo 

/•*/2 „2 

< C7i(e2 + /i(t)2)(l + t)T+V2 



2 rt 



l + t)~i/ e"Ards + e"2M / (i + t_s)-ids 

Jt/2 



which completes the proof of the lemma. 



□ 
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